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I. INTRODUCTION 
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body system. This is either due to high precision that calls for an inclusion of all possible 



Interpretation of many modern experiments requires a covariant description of the two- 



corrections to a standard (possibly nonrelativistic) approach or due to the high energies and 
momenta involved in the processes investigated. 

In the spirit of a local quantum field theory the starting point of a relativistic covariant 



description of bound states of two particles is the Bethe-Salpeter (BS) equation. However, 



despite the obvious simplicity of two-body systems, the procedure of solving the BS equation 
£T) ', encounters difficulties. These are related to singularities and branch points (cuts) of the 

in ■ 

■^J- ■ amplitude along the real axis of the relative energy in Minkowski space. Therefore, up to now 



the BS equation including realistic interaction kernels has been solved either in Euclidean 



OO 

space within the ladder approximation or utilizing additional approximations of the equation 
itself. 

Unfortunately, our understanding of the mathematical properties of bound states within a 
relativistic approach is far from being perfect. In mathematical terms the BS equation itself 
is a quite complicated object, and the technical problem of solving it is still a fundamental 
issue. Consequently there are very few successful examples of solving the BS equation for 
fermions including realistic interactions. For example, in Ref. 1J the BS equation for spinor 
particles was regularly treated by using a two-dimensional Gaussian mesh. That series of 
studies revealed a high feasibility of the BS approach to describe nucleon-nucleon interactions, 
in particular, processes involving the deuteron. However, it should be mentioned that the 
algorithms exploiting the two-dimensional meshes are rather cumbersome and require large 
computer resources. In addition, the numerical solution is obtained as two-dimensional arrays 
which are quite awkward in practice when computing matrix elements and in attempts to 
establish reliable parameterizations and possible analytical continuations of the solution back 
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to Minkowski space. Therefore, it is necessary to provide a method for solving BS equations 
that would feature a smaller degree of arbitrariness. 

In the present paper we suggest an efficient method to solve the BS equation for fermions 
involving interaction kernels of one-boson exchange type supplemented by corresponding form 
factors. It is based on hyperspherical harmonics used to expand partial amplitudes and 
kernels. We show that this novel technique provides many insights into the BS approach. 
The current study is partially stimulated by the results reported in Ref. 2j. We explore the 
structure of X S and 3 Si — 3 D\ bound states for different couplings and study the details of 
the convergence of solutions and corresponding eigenvalues. In particular, on the basis of 
the proposed method for solving the BS equation it becomes possible to analyze the specifics 
of the problem related to the stability of bound states in the BS approach. Besides, the 
hyperspherical expansion provides an effective parameterization of the amplitude, which is 
extremely useful in practical calculations of observables and in theoretical investigations of the 



separability of the BS kerne 
of method may be found in 



with one-boson exchange interaction. The detailed description 
3,. 



II. OVERVIEW OF THE METHOD 

The BS equation for two spinor particles interacting via one-boson-exchange potentials for 
vertex Q{p) being a 4 x 4 matrix in spinor space has the form 

0(p) = J j0p v( P , k) r(i) s(h) g(k) s(k 2 ) f(2), (i) 

where the propagator V(p, k) for scalar and pseudoscalar exchange mesons is 

V (P> k ) = 7 2^ ~ ( 2 ) 

(p — k) 2 — n l + %e 
while for spinor particles the corresponding propagators are 

k\+ra ~ tl N _ ^r,,, , T k 2 — m 



s ( k J = >2 9 , • . s (^) = cs(k 2 y c 



k\ — m 2 + is ' " ~ k\ — m 2 + is ' 

where C = i r f° r j 2 is the charge conjugation matrix, k± t 2 and px,2 are the four-momenta of 
the constituent particles, k and p are the corresponding relative momenta. The meson ver- 
tices T(l,2) are determined by the effective interaction Lagrangians urged to describe the 
considered fermion system. For a two nucleon system, within the one-boson exchange ap- 
proximation, these vertices are T(l) = 1; T(2) = —1 for scalar and T(l) = 75; T(2) = —75 
for pseudoscalar couplings, respectively. Each interaction vertex T is augmented with cut-off 
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monopole form factors F(q 2 ) = A 2 / (A 2 — q 2 ) where A are free parameters. Note that the 
coupling constant g in eq. §I§ is purely imaginary for the pseudoscalar mesons else purely 
real. 

The first step is the expansion of the vertex function Q{po,p) into spin-angular harmonics 

G(po,p) = ^2g*(poAp\)Z*(p)- (3) 

a 

For specific bound states with given quantum numbers only some basis matrices contribute 
to the vertex function G(p) . E.g., for the 1 5'o state only four matrices are relevant to describe 
the amplitude, while in the 3 Si — 3 D\ channel eight basis matrices are needed. In the 1 Sq 
channel the basis is 

Up) = ^75; Up) = ^=7075; Up) = -^= {j ^ioi 5 ; Up) = ~4^^b. (4) 

V167T V167T ylQit \P\ V167T \P\ 

By employing the Pauli principle and the charge conjugation operation one obtains [3j that 
in the 1 5'o channel the component g<± is of the odd parity while the remaining g%, g^ are 
of the even parity. As mentioned, the BS amplitude is a mathematically complicated object 
and it is more convenient to considered it in Euclidean space, where the analytical properties 
of the amplitude become simpler and more transparent. For convenience, in Euclidean space 
we redefine the odd partial components g^ for the 1 5'o channel as g^ — > ig^. Then the Wick 
rotated BS equation Q reads 

9n(P4,\p\) = g 2 J dtt p J ^-^£(£4, |fc[) ^ _ + ^ ^A nm (p, k)g m (h, k), (5) 

where m,n = 1 ... 4 for the 1 Sq and the scalar part S(k^, \ k\) of the two spinor propagators 
is defined as 

S(h,\k\) = — ^ . (6) 

(k 2 + m 2 -^) 2 + M 2 k 2 

Next step is the expanding of the interaction kernel into hyperspherical harmonics 4] 
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(p — k) 2 + fi 2 t-* 1 n + 



2tt 2 Y] —^—V n (p, k)Z n i m (xp, P , (p P )Z* lm (xk, k , <f>k), (7) 
/ — ' n + 1 

nlm 



A, - A_ 



with 



'2 2 ^ (n + l)(n - l)U\ 2 ^ , 



Z n im(x,0,(j)) = X nl (x)Y lm (6,(j)); X nl (x) = \ ? — ——-—sin xC n +;(cosx), 

7r [n + I + 1 ) ! 



4 



where Yi m (9, <p) are the familiar spherical harmonics, and C l ^\ are the Gegenbauer polynomi- 
als. The resulting system of equations reads 

9n(p4, \p\) =g 2 J k dkS ^* kdXk S{k 4 , \k\)Wi n (p,k, 

XpiXk) / j Q"n,m 

(k 4 ,k)g m (k 4 ,k), (8) 



where 



Wi n (p,k,Xp,Xk) = ^2f^l v i(Pi ~ k ) x u n (Xp)Xii n {Xk)- 



+ 

Furthermore, the partial vertex functions g n are expanded over the basis X n i(x P ) as 

oo oo 

01,2 (P4, \P\) = ^2gi2(P) X 2j-2,o(X P ); 03(P4, \p\) = dip) X 2j -X,l(Xp)] (9) 
3=1 3=1 

oo 

9i(p4, \p\) = ^2gi(P)x 2j ,i(xp)- (io) 

3=1 

Placing eqs. (Q and (fTUj) in to eq. and performing integration over hyper-angles ana- 
lytically, the initial system of four-dimensional integration equations reduces to a (infinite) 
system of only one-dimensional equations. Eventually, by a proper choice of numerical inte- 
gration method (the Gaussian formula, in our case) one easily transforms the latter system of 
integral equations in to an ordinary algebraic system of linear equations. Obviously, further 
solving procedure is straightforward (see for details Ref. [sj]). 



III. RESULTS 



By the above procedure we solved the BS equation for a system of two spinors with equal 
masses in 1 Sq state interacting via exchanges of scalar and pseudoscalar mesons. Herebelow, 
for the sake of brevity, we present results only for scalar exchange mesons; more detailed 
results can be found in [3|. Note that the accuracy of the proposed method depends up 
on the dimension Nq of the Gaussian mesh used in calculations and on number of terms 
M max used in the decomposition (jUJ) and ffTUj) . Also the magnitude of the coupling constants 
governs the existence of the solution itself; at some critical values of the coupling constants 
the solution of the BS equation may not exist at all. 

As an illustration of the stability of the numerical procedure, in Table [I] we present results 
for the masses of the bound state M(g 2 ) depending on the Gaussian mesh Nq and M max . 
Calculations have been performed for the 1 So state, with a scalar meson exchange of mass \i 
for two values // = 0.15 GeV/c 2 and /x = 0.5 GeV/c 2 ; the constituent particles (nucleons) have 
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been taken with equal masses m = 1.0 GeV/c 2 for simplicity. Results presented in Tabled 
clearly demonstrate that the approximate solution converges rather rapidly, and already at 
M max ~ 4 — 5 and Ng = 64 the method provides a good solution of the system. 



<? 2 = 15 


// = 0.15 GeV/c 2 


At = 0.5 GeV/c 2 


-'"max 


N G = 32 


Ng = 64 


N G = 96 


N G = 32 


N G = 64 


N G = 96 


1 


1.9399 


1.9399 


1.9399 


1.9984 


1.9984 


1.9984 


2 


1.9370 


1.9370 


1.9370 


1.9982 


1.9982 


1.9982 


3 


1.9368 


1.9368 


1.9368 


1.9982 


1.9982 


1.9982 


4 


1.9368 


1.9368 


1.9368 


1.9982 


1.9982 


1.9982 




<? 2 = 30 


// = 0.15 GeF/c 2 


At = 0.5 GeV/c 2 


-'"max 


Ng = 32 


N G = 64 


N G = 96 


N G = 32 


N G = 64 


iV G = 96 


1 


1.7932 


1.7910 


1.7905 


1.9167 


1.9142 


1.9137 


2 


1.7897 


1.7875 


1.7871 


1.9152 


1.9127 


1.9122 


3 


1.7896 


1.7874 


1.7870 


1.9152 


1.9127 


1.9122 


4 


1.7896 


1.7874 


1.7870 


1.9152 


1.9127 


1.9122 



TABLE I: Dependence of the bound state masses on the M max and on the Gaussian mesh used in 
actual numerical calculations. 

The solution of the BS equation, eq. ([3]), can be obtained from the known numerical 
values of the partial components <fc. In Fig. [TJ we present the behavior of the partial co- 
efficients in the 1 So channel g{(p), j = 1...4, eq. (jUJ), as a function of the euclidian 
relative momentum p = yp\ + p 2 . Calculations have been performed for a bound sys- 
tem with M = 1.937 GeV/c 2 . The actual parameters used in numerical calculations are 
Ng = 96, M max = 4 and g 2 = 15. Closed squares correspond to g\, closed circles - g\ 
multiplied by 10, triangles -g\ multiplied by 100, open circles -g\ multiplied by 1000. It 
can be seen that at large p each function decreases as inverse powers of p, which allows 
for a relatively simple parametrization of the result; the solid lines correspond to such a 
fit. An analysis of the obtained results shows that the numerical solution is rather sensi- 
tive to the magnitude of the coupling constants. Moreover, it has been found that there 
are some critical values of the coupling constants, g 2 v , above which the solution of the BS 
equation does not exist (in absence of the cut-off form factors). At values of the coupling 
constants close the their critical values g 2 ~ g 2 r the numerical solution becomes unstable and 
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FIG. 1: Left panel: The coefficients g\, j = 1 . . . 4, eq. ([3p, as functions of the euclidian momentum 
P 2 = p\ + P 2 - Right panel: The mass of the bound state 1 Sq as a function of the cut-off parameter 
A at different values of the coupling constant g. At large g > g cr the solution is strongly dependent 
on the cut-off parameter, i.e. becomes rather unstable. 

strongly dependent on the Gaussian mesh, M max and A. Such a situation is illustrated in 

Fig. [H right panel, where we present the dependence of 
the mass M on the values of the cut-off parameter A at 
different coupling constants g 2 . It is clearly seen that at 
relatively low values of g 2 below its critical value (g 2 r ~ 40) 
the solution is practically independent on A, i.e. it con- 
verges rather rapidly. Contrarily, at g 2 > g 2 T the solution 
M(g 2 ) becomes manifestly dependent on A, i.e., it can dis- 

"'"5 10 15 20 25 30 appear at all. Such a behavior of the solution at g 2 ~ g 2 r 
2 

9 exactly reproduces the peculiarities of the well-known col- 

lapse phenomenon for potentials like —a/r 2 in nonrela- 
FIG. 2: The mass of the bound tivistic quantum mechanics. The magnitude of the crit- 
state M in the 1 Sq channel within ical value g 2 r may be estimated analytically and it turns 
different models, non relativistic out to be g 2 v ~ Arc 2 (see Ref As mentioned, the 

(NR), light front (LF) dynamics and obtained solution has a rather simple dependence on the 
present (BS) approach. euclidian momentum p and allows for a simple analytical 

parametrization of the solution, which can be exploited further to perform a Wick rotation 
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back to Minkowsky space. In this context is is instructive to compare our solution with ones 
obtained within other known approaches. In Fig. [2] a comparison of our solution (solid line) 
with the results of non relativistic calculations (dotted line) as well as with the ones within the 
Light Front dynamics (LF) [2] (dashed line) is presented. As expected, at low values of the 
coupling constant different approaches provide similar results. As g increases the difference 
becomes more significant, reflecting the role of relativistic effects. The difference between LF 
and BS approaches is due to different treatment of the vertex function, namely within the LF 
the vertex function consists on three components, while within the BS formalism four par- 
tial components describe the solution. Obviously, the role of the fourth component increases 
with g 2 . 



IV. CONCLUSION 



We generalize a method based on hyperspherical harmonics to solve the homogeneous 
spinor-spinor Bethe-Salpeter equation in Euclidean space. To do so, we introduce a new basis 
of spin-angular harmonics, suitable to expand the Bethe-Salpeter vertex into four-dimensional 
hyperspherical harmonics. We obtain an explicit form of the corresponding system of one- 
dimensional integral equations for the partial components and formulate a proper numerical 
algorithm to solve this system of equations. The BS vertex functions are studied in detail for 
the 1 S'o and 3 Si — 3 D% bound states with scalar, pseudoscalar and vector meson exchanges. 
Our results are in a good agreement with calculations within the non relativistic and Light 
Front Dynamics approaches. Within the novel method the effectiveness of the numerical 
procedure is analyzed for the scalar, pseudoscalar and vector meson exchanges and conditions 
for stability of the solution are established. 

An advantage of the method is the possibility to present the numerical solution in a reliable 
and simple analytical parameterized form, extremely convenient in practical calculations of 
matrix elements within the BS formalism and for analytical continuation of the solution back 
to Minkowski space. The method allows for a covariant description of two-body systems, such 
as the deuteron, positronium and the variety of known mesons, as relativistic bound states 
and, by solving the inhomogeneous BS equation, to describe the scattering states. 
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